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1  Introduction 


Purpose 

Wave  measttrementa  are  required  to  observe  and  ultimately  predict  tbe 
response  of  tbe  ocean  to  dynamic  loading.  While  remote  sensing  tech¬ 
niques  may  eventually  inovide  continuous  spatial  definition  of  wave  en¬ 
ergy  for  large  tegi<ms.  current  engineering  practice  relies  on  point 
measurements  firom  wave  gauges.  The  U.S.  Army  Corps  of  Engineers’ 
Field  Wave  Oa^ng  Program  (PWQP)  operates  a  rutionwide  network  of 
wave  gauges  and  serves  as  the  Corps*  central  access  point  for  all  U.S. 
wave  data.  The  program  also  supports  a  variety  of  activities  to  enhance 
die  quality  and  efficiency  of  wave  data  collecdtm,  processing,  and  distribu- 
titm.  This  is  the  diird  in  a  series  of  FWOP  reports  descriUng,  through 
case  studies,  the  process  of  selecting  sites  for  wave  gauges. 

The  selecdcm  of  sites  for  measuring  wave  energy  in  a  particular  envi¬ 
ronment  is  oftmi  a  heuristic  exercise  wherein  engineering  judgement  is 
used  to  balance  tbe  sometimes  conflicting  requirements  of  defirditg  the 
processes  (physics),  ensudng  data  capture  Qogistics),  and  managing  re¬ 
sources  (economics).  Balancing  these  aspects  may  never  yield  to  purely 
arudytic  soludon,  but  tools  can  be  devel(^>ed  to  quantify  befits  and 
costs.  In  O’Reilly  and  McOehee  (in  preparadon)  a  ted^que  borrowed 
firom  metallurgy-simulated  annealing  is  described  that  optimizes  the  infor¬ 
mation  ccmtent  of  a  network  within  a  region  (assumed  to  have  uniform  in¬ 
cident  conditions  at  its  outer  boundary).  Basco  and  McOehee  (1990)  walk 
through  all  tbe  steps  involved  in  selecting  the  exact  gauge  location  once  a 
general  site  has  been  identified.  This  pqier  illustrates  use  of  two  different 
wave  transformation  models  to  optimize  gauge  sites  in  a  local  gauge  net¬ 
work  for  a  harbor  oscillation  study.  The  harbor  under  consideration  is 
Kahului  Harbor,  Hawaii. 


Background 

As  with  the  other  Hawaiian  Islands,  Maui  is  almost  totally  dependent 
on  ocean  transportation  for  its  basic  supplies,  and  Kahului  Harbor  is  the 
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only  commercial  deep-draft  harbor  on  this  island.  In  order  to  handle  the 
growing  cargo  tonnage  generated  by  the  growth  in  population  and  expand¬ 
ing  economy,  and  to  keep  pace  with  the  technological  and  operational 
changes  occurring  in  the  maritime  industry,  modifications  to  the  present 
port  facilities  are  being  planned  (Hawaii  Department  of  Transportation 
(HDOT)  1989).  HDOT  requested  that  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  Coastal  Engineering  Research  Center  (CERC)  conduct 
a  study  of  the  existing  conditions  at  Kahului  Harbor  and  make  recommen¬ 
dations  to  maximize  the  use  of  available  areas  for  harbor  expansion  and 
navigational  improvements. 

Although  existing  swell  conditions  at  Kahului  Harbor  have  not  resulted 
in  complete  closure  of  the  harbor,  they  have  hampered  the  efficient  use  of 
some  berths.  Though  unlikely,  there  is  concern  that  supply  lines  to  the  is¬ 
land  might  be  cut  off  if  the  entrance  channel  were  temporarily  blocked 
(HDOT  1989).  It  is  not  clear  whether  the  existing  problems  are  due  to 
swell  entering  the  harbor  or  lower  frequency  ‘seiche*  motions,  though  sig- 
nifrcant  swell  can  be  observed  throughout  the  harbor  when  northerly 
storms  occur.  Operational  problems  in  many  harbors  are  due  to  seiche 
(also  called  surge  or  resonance)  which  is  characterized  by  a  concentration 
of  energy  in  standing  wave  modes  at  the  resonant  frequencies,  typically  be¬ 
tween  0.001  and  0.01  Hz  for  small  harbors.  Energy  at  these  frequencies  is 
amplified  within  the  harbor,  relative  to  outside  the  harbor,  whereas  harbor 
motions  at  other  frequencies  are  suppressed.  In  harbors  the  size  of 
Kahului,  the  lowest  frequency  resonances  may  be  driven  by  atmospheri¬ 
cally  excited  internal  and  shelf  waves,  whereas  higher  frequency  modes 
are  usually  driven  by  nonlinear  interactions  between  swell  and  sea.  The 
forcing  and  damping  of  resonant  modes  may  differ  from  harbor  to  harbor 
depending  on  the  harbor  geometry,  surface  wave  climatology  and  general 
oceanographic  properties,  (i.e.,  stratification  and  wind  field)  in  the  sur¬ 
rounding  offshore  area.  It  is  not  presently  possible  to  quantitatively  pre¬ 
dict  the  magnitude  of  harbor  seiches  either  numerically  or  with  physical 
models,  owing  to  the  complexity  of  the  seiche  excitation  and  damping  pro¬ 
cesses.  In  situ  observations  are  essential  to  the  characterization  of  barter 
seiche  and  to  provide  the  unknown  values  of  parameters  necessary  for 
modeling  studies. 


Kahului  Harbor  Site  Description 

Kahului  Harter  is  situated  at  the  apex  of  a  large  V-shaped  bay  formed 
by  an  indentation  on  the  north  shore  of  the  island  (Figure  1).  Offshore  ba¬ 
thymetry  at  this  site  is  complicated  (Figure  2).  East  of  the  barter  the  ba¬ 
thymetry  is  characterized  by  shallow  areas  (<  3  m)  which  extend  nearly 
1-2  km  offshore  whereas  to  the  west  the  bottom  slope  is  steeper  and  rela¬ 
tively  constant.  Directly  offshore  of  the  entrance  channel  the  depth  is 
about  10  m. 
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Figure  1.  Kahului  Harbor  Site,  Maui.  Hawaii 

Overall  harbor  dimensions  are  approximately  1  km  by  1  km,  although 
depths  along  the  unused  western  and  southern  perimeter  are  only  a  few 
meters  (Figure  3).  The  present  operational  harbor  basin  defined  by  the  10* 
to  11-m  depth  contour  (Figure  3)  is  approximately  730  m  long  and  600  m 
wide.  The  -200-m-wide  entrance  chaimel  faces  northward  and  is  bounded 
by  two  rubble-mound  breakwaters.  Although  under  most  wave  conditions 
these  breakwaters  successfully  protect  the  harbor,  energetic  nordierly 
swell  teaches  all  piers  and  sometimes  interferes  with  harbor  operations 
(HDOT  1989). 
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Figure  2.  Offshore  and  harbor  bathymetry  at  Kahului,  Maui.  Dashed  line  indicates  the 
near-field  HARBD  numerical  model  grid  boundary  for  Test  1  and  the  dotted 
line  is  the  boundary  for  ail  other  tests.  0  (top  of  figure)  is  the  incident  wave 
direction.  The  site  selected  for  the  offshore  slope  array  is  shown  by  a 
crossed  square 
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Figure  3.  Harbor  bathymetry  and  sensor  locations.  The  solid  symbols 
indicate  positions  for  which  HARBO  model  results  are  shown 
in  subsequent  figures  and  the  lettered  circles  are  the  harbor 
sites  selected  for  field  data  collection 
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Approach 


General 

The  gauging  plan  arrived  at  in  this  report  is  one  element  of  a  broader 
investigation  planned  for  Kahului  Harbor. 

CBRC  has  proposed  a  comprehensive  study  to  address  the  effect  of  var¬ 
ious  modifications  to  the  harbor  on  amplitude  and  incidence  of  occurrence 
of  both  swell  and  harbor  seiche  throughout  the  harbor.  The  approach  is  to 
utilize  a  three-dimensional  physical  model  and  a  numerical  harbor  oscilla¬ 
tion  model  (the  HARBD  model)  to  simulate  tiie  response  of  the  existing 
harbor  to  gravi^  and  infregravity  energy.  These  models  will  be  calibrated 
and  verified  with  field  measurements.  A  directional  array  placed  outside 
the  harbor  will  measure  incident  waves  that  excite  the  harbor.  Four  non- 
directiorud  gauges  placed  throughout  the  harbor  wUl  measure  the  response 
of  the  harbor.  Once  the  models  have  been  validated,  alternative  plans  will 
be  compared. 

Two  additional  models  will  be  utilized  to  define  the  statistical  wave  cli¬ 
mate  outside  the  harbor  so  that  the  probability  of  occurrence  of  different 
forcing  conditions  can  be  estimated.  The  first  is  WISWAVE,  the  model 
used  by  the  CERC  Wave  Information  Study  (WIS)  (Corson  et  al.  1986)  -  a 
hindcast  of  deep-ocean  wave  conditions  at  discrete  stations  based  on  long¬ 
term  meteorological  observations.  Completed  for  the  20-year  interval 
from  1956-1975,  the  WIS  hindcast  is  scheduled  in  1996  for  updating 
through  199^  which  will  cover  the  period  when  the  Kahului  Harbor 
gauges  are  :^«rating.  The  updated  hindcast  will  be  validated  against  a  di¬ 
rectional  wave  gauge  located  approximately  40  miles  northwest  of  Maui. 
This  gauge  was  deployed  in  Jamiary  1993  by  the  National  Data  Buoy  Cen¬ 
ter  and  was  funded  by  the  FWOP  as  part  of  a  joint  effort  between  the 
Corps  and  the  State  of  Hawaii  to  define  the  wave  climate  for  the  islands. 

The  final  model  requited  is  a  site-specific  wave  transformation  model 
that  will  predict  wave  conditions  outside  the  harbor  using  the  validated 
WIS  hindcast  as  input  Due  lo  the  comidex  bathymetry  approaching  the 
harbor,  this  model  will  incorporate  both  refraction  and  difihection  effects. 
The  wave  transformation  model’s  ouqmt  will  be  validated  through 


6 


ChaptarZ  Approach 


comparisons  with  the  directional  array  used  for  input  to  the  harbor  oscilla¬ 
tion  model.  Figure  4  illustrates,  and  Table  1  summarizes,  the  progression 
of  models  to  be  used  to  provide  the  statistical  predications  of  the  harbor 
response. 

To  design  the  Kahului  gauge  network,  two  numerical  models  were  run 
in  advance  of  the  gauge  placement  to  indicate  the  optimum  locations  for 
the  gauges  -  the  HARBD  harbor  oscillation  model  (Chen  1986)  and  the  re¬ 
fraction/diffraction  (R/D)  wave  transformation  model  developed  at  the 
Scripps  Institution  of  Oceanography  (SIO)  (O’Reilly  and  Guza  1993). 
Output  from  the  models  was  examined  to  support  or  reject  sites  selected 
through  the  more  traditional  “inspection”  approach. 

While  quantitative  results  from  the  uncalibrated  numerical  models  are 
not  always  reliable,  valuable  qualitative  information  can  be  derived  by 
running  the  models  prior  to  obtaining  the  measurements.  The  objectives 
were  to  predict;  (a)  probable  modes  of  harbor  oscillation;  (b)  sensitivity 
of  the  harbor  system  to  variations  in  geometry  and  bathymetry;  and  (c)  lo¬ 
cations  outside  the  harbor  that  are  strongly  influenced  by  local  bathyme¬ 
try,  i.e.,  where  conditions  would  be  unrepresentative  of  the  incident 
conditions  for  the  harbor  models. 

The  principal  concern  expressed  by  the  HDOT,  and  thus  the  emphasis 
for  the  gauge  selection  process,  was  the  potential  response  of  conunercial 
vessels  to  harbor  oscillations  while  moored  at  their  berths.  Thus,  while 
HARBD  is  capable  of  modeling  infragravity  and  gravity  waves,  it  was  run 
here  to  identify  the  harbor  seiche  modes.  It  is  proposed  to  also  be  cali¬ 
brated  for  swell  during  the  comprehensive  study.  The  R/D  model  was  run 
at  a  single  representative  gravity  (0.058  Hz,  or  17  sec)  wave  approaching 
from  two  directions  to  investigate  the  variability  of  swell  conditions  out¬ 
side  the  harbor. 


Physics 

In  order  to  run  the  models  in  advance,  the  best  available  information  on 
bathymetry  inside  and  outside  the  harbor  was  utilized.  Additional  experi¬ 
ence  with  the  models,  and  examination  of  the  measurements,  will  help  de¬ 
termine  if  additional  or  more  updated  bathymetric  data  are  required. 

Other  assumptions  were  made  about  the  reflectivity  of  harbor  boundaries 
and  the  configuration  of  fill  underneath  piers  that  will  have  to  be  verified. 

While  all  of  the  numerical  models  above  are  spectral,  they  are  linear  in 
frequency;  i.e.,  they  can  only  predict  changes  in  energy  or  direction 
within  each  single  frequency  band  of  the  spectrum.  Observations  in  na¬ 
ture  indicate  that  nonlinear  energy  transformation  from  gravity  waves  to 
infragravity  waves  is  responsible  for  excitation  of  harbor  seiche  during 
storm  conditions.  The  proposed  physical  model  has  the  potential  to  dupli¬ 
cate  this  process,  though  questions  remain  on  the  scale  effects.  As  part  of 
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Rgure  4.  Progression  of  numerical  models  proposed.  Full  ocean  scale  WIS  hindcast  output  drives  wave  transformation  model, 
put  of  wave  transformation  model  is  input  for  HARBO  and  3-D  physical  models 


table  1 

Summary  of  Modala 


Modal 

Seale 

Grid 

Input 

Output 

Validate 

WIS  numarieal 
wave  generated 
model 

North  Pacific 

0  100  km' 

Atmospheric 
pressure,  20  -f 
years 

2-D  spectra@ 
Phase  1  stations 

NDBC  buoys 

R/0  numerical 

wave 

transformation 

model 

Deep  ocean  to 

harbor 

entrance 

0  100  m' 

WIS  Phase  1; 
and/or  NDBC  buoy 

2-D  spectra® 
harbor  entrance 

Directional 
array  outside 
harbor 

HARBD  numerical 
harbor  oscillation 
model 

Harbor  and 
vicinity 

O  10  m' 

R/D  model  output; 
and/or  directional 
anay 

Amplification 
spectra  inside 
harbor 

Nondirectional 
gauges  in 
harbor 

3-D  physical 
harbor  model 

Harbor  and 
vicinity 

Continuous 

R/D  model  output; 
and/or  directional 
array 

2-D  spectra 
inside  harbor 

Nondirectional 
gauges  in 
harbor 

'  0 10  m  is  shorthand  tor  *on  the  order  of  10  m.' 


the  study,  empirical  relations  between  swell  energy  and  seiche  energy  will 
be  obtained  from  the  wave  measurements.  This  relationship  will  be  uti¬ 
lized  to  verify  similarity  in  the  physical  model  and  to  develop  the  statisti¬ 
cal  predictions  of  harbor  response. 

The  HARBD  model  assumes  that  all  energy  is  either  reflected  or  ab¬ 
sorbed  at  boundaries,  but  no  transmission  occurs  through  boundaries.  A 
method  of  representing  the  permeable  walls  will  be  described,  but  no  tests 
of  the  accuracy  of  this  ‘‘simulated  permeability”  have  been  conducted. 
Physical  models,  properly  scaled  and  calibrated,  will  also  be  used  to  simu¬ 
late  the  transmission  of  energy  through  and  over  the  structure. 


Logistics 

The  models  described  above  require  measurements  of  wave  energy  as  a 
function  of  frequency  (1-D  spectrum)  or  frequency  and  direction  (2-D 
spectrum).  Two  generic  types  of  wave  gauges,  buoys  and  bottom- 
mounted  pressure  gauges,  are  generally  utilized  to  obtain  these  data 
(Hemsley,  McGehee,  and  Kucharski  1^1).  A  surface  following  buoy  was 
selected  for  the  deep-ocean  WIS  calibration  site,  and  bottom-mounted 
pressure  gauges  for  the  harbor  network.  Single  pressure  gauges  are  used 
to  obtain  1-D,  or  nondirectional  spectra,  and  2-D  spectra  are  obtained 
from  an  array  of  four  sensors  on  a  single  frame. 

Bottom-mounted  instruments  can  be  either  self-contained  or  cable  con¬ 
nected  to  shore.  Armored  cables  provide  power  to  the  underwater  sensor 
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and  carry  the  signal  to  electronic  systems  on  shore,  where  they  can  be 
stored  and/or  telemetered  (via  radio  or  telephone)  to  a  central  collection 
computer.  A  real-time  reporting,  cable-connected  network  was  selected  for 
the  harbor  gauges  to  improve  chances  of  data  capture,  and  to  provide 
HDOT  with  the  potential  for  developing  a  wave  monitoring  system  for  use 
in  harbor  operations. 

System  acquisition  and  installation  costs  are  affected  by  the  length  of 
the  cable;  therefore,  it  is  desirable  to  minimize  the  distance  from  the  shore 
electronics  station  to  the  sensor.  Practical  locations  for  the  shore  station 
must  have  protection  from  elements  and  tampering,  and  access  to  1 10-V 
electrical  power  and  telephone  lines.  Locations  were  identified  during  pre¬ 
vious  site  visits  at  the  eastern  comer  of  the  harbor  and  the  outer  end  of 
pier  2  (Figure  2).  An  upper  limit  on  the  cable  length  is  approximately 
2  km  before  expensive  in-line  ampliHers  are  required,  so  all  sensor  sites 
must  be  within  2  km  of  either  shore  station  site. 

While  the  underwater  sensors  are  subject  to  electronic  failure  and 
water  leakage  problems,  experience  has  shown  ship  activity  to  be  the  prin¬ 
cipal  hazard.  Activities  around  the  harbor  that  threaten  either  the  gauges 
or  their  cables  are  anchoring  and  dragging  of  barge  tow  bridles.  For  this 
reason,  instruments  and  cable  routes  should  avoid  anchorages  and  the  prin¬ 
cipal  navigation  channels  -  particularly  near  the  entrance  where  tow  bri¬ 
dles  are  adjusted. 

Installation  of  the  bottom-mounted  sensors  requires  divers.  As  installa¬ 
tion  depths  increase,  significant  time  penalties  on  diver  bottom  time  accu¬ 
mulate.  Another  depth  limitation  is  the  attenuated  frequency  response  of 
pressure  sensors  with  increasing  depth.  Both  factors  dictated  a  maximum 
sensor  depth  of  about  20  m.  For  additional  information  on  gauges  and 
gauge  site  selection,  see  Hemsley,  McGehee,  and  Kucharski  (1991). 

A  major  factor  influencing  the  site  selection  process  is  the  areas  of  in¬ 
terest  from  the  perspective  of  the  port  activities.  Current  operations  are 
principally  in  the  eastern  comer  and  northeastern  side  of  the  harbor.  Fu¬ 
ture  expansion  options  (2010  Master  Plan,  HDOT  (1989))  are  the  southern 
side,  west  of  Pier  2,  and  the  northwestern  side  of  the  harbor.  Proper  vali¬ 
dation  of  the  models  for  these  areas  in  their  current  configuration  is  im¬ 
portant  for  establishing  confidence  in  predictions  for  new  configurations. 
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3  HARBD  Numerical  Model 


HARBD  Model  Description 

The  sometimes  complex  spatial  stracture  of  harbor  seiches,  with  min¬ 
ima  and  maxima  of  sea  surface  elevation  (i.e,  nodes  and  antinodes),  is  de¬ 
termined  by  the  harbor  geometry  and  bathymetry.  The  spatial  structure  of 
the  resonant  modes,  within  the  frequency  range  of  likely  seiches  at 
Kahului  Harbor  (0.001  -  0.01  Hz)  was  obtained  numerically  using  the 
Corps  HARBD  model  (Chen  1986,  Chen  and  Houston  1987).  The  numeri¬ 
cal  model  is  the  solution  of  a  boundary  value  problem  for  linear  gravity 
waves  in  a  domain  of  arbitrary  shape  and  variable  depth  (Chen  and  Mei 
1974,  Houston  1981).  The  model  domain  is  divided  into  a  variable-depth 
near  field  and  a  constant-depth  far  field,  separated  by  an  artificial  semicir¬ 
cular  boundaiy  where  solutions  are  matched.  The  solution  in  the  near 
field  is  found  using  finite  element  techniques  and  the  far-field  solution  is 
analytic. 

Dissipative  effects  of  bottom  friction  and  boundary  absorption  can  be 
included  in  the  model,  and  the  sensitivity  of  the  modal  spatial  structure  to 
these  sources  of  energy  losses  has  been  investigated.  Bottom  friction  is  as¬ 
sumed  to  be  linearly  proportional  to  the  horizontal  velocity  and  is  pre¬ 
scribed  in  terms  of  a  spatially  variable  bottom  friction  coefficient  fi  where 
the  inviscid  solution  is  obtained  using  everywhere.  Side-wall  bound¬ 
ary  absorption  is  prescribed  in  terms  of  a  reflection  coefficient  R,  where 
Rsl  implies  perfect  reflection.  HARBD  does  not  allow  wave  transmis¬ 
sion  at  boundaries,  so  if  waves  impinging  on  the  breakwaters  are  only  par¬ 
tially  reflected,  the  remainder  is  absorbed. 

Undesirable  wave  reflection,  refraction,  and  scattering  occur  at  the 
depth  discontinuity  along  the  semicircular  boundary  between  the  variable- 
depth  near  field  and  the  constant-ctepth  far  field.  Effects  of  these  depth 
discontinuities  are  reduced  by  placing  the  semicircular  boundary  where 
the  actual  bathymetry  is  relatively  snoooth  and  planar  and  selecting  a  con¬ 
stant  far-field  depth  matching  the  measured  bathymetry  over  a  large  por¬ 
tion  of  the  boundary  (Okihiro,  Guza,  and  Seymour,  in  preparation). 

Nonzero  bottom  fiiction  can  also  be  specified  near  the  boundary  between 
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the  coast  and  the  semicircular  boundary  to  further  reduce  the  bogus  ef¬ 
fects  of  these  depth  discontinuities. 


Finite  Eiement  Grids 


The  bathymetric  grid  used  in  the  near-field  model  solution  was  con¬ 
structed  from  bathymetric  surveys  obtained  from  the  Corps.  HDOT,  and 
the  National  Ocean  Service/National  Oceanic  and  Atmospheric  Adminis¬ 
tration.  Very  little  data  was  available  for  the  shallow  areas  (i.e.,  <  2  m) 
both  inside  and  outside  the  harbor.  The  sensitivity  of  the  model  results  to 
the  specification  of  depths  in  these  shallow  areas  and  the  magnitude  of  the 
discrepancies  between  the  model  results  and  observation  can  be  used  to  de¬ 
termine  whether  additional  bathymetry  data  should  be  incorporated  into 
the  model  grids  for  future  work. 

Mean  sea  level  was  used  in  all  the  model  tests.  The  constant  depth  in 
the  far  field  was  chosen  to  be  14  m.  The  triangular  elements  had  diroen- 
sitms  of  15-30  m  (less  than  1/6  of  the  local  wavelength).  Model  tests 
were  run  on  the  C^ce  of  Naval  Research  Cray  Y-MP  and  the  Center  for 
Coastal  Studies  (at  SIO)  Sun  670MP. 
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4  HARBD  Model  Tests  and 
Results 


HARBD  Modal  Tests 

The  HARBD  model  tests  are  summarized  in  Table  2.  Owing  to  limita- 
dam  on  available  cmnpoter  resources,  model  results  were  obtained  for  a 
single  incident  wave  direction,  approximately  18  deg  east  of  north  (Fig¬ 
ure  2)  and  for  frequencies  less  than  0.012  Hz.  Limited  model  tests  with 
varying  incidmit  wave  directions  showed  that  while  die  magnitude  of  the 
resonant  peaks  may  vary  somewhat  as  a  foncdon  of  incident  wave  direc- 
dtm,  die  resonant  firequendes  are  not  highly  smisidve  to  direcdon.  Based 
on  previous  work  with  a  similar  sized  hartxir  (Okihiro,  Ouza,  and  Sey¬ 
mour,  in  preparadon)  most  of  die  seiche  oiergy  is  expected  to  be  at  fre¬ 
quencies  less  than  ~0.012  Hz. 


Tibte2 
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Inviscid  Tests 


The  first  set  of  model  tests  (Tests  1  end  2,  Tsble  2)  neglected  bottom 
friction  (i.e.,  baO)  and  assumed  complete  reflection  (Ral)  along  the 
coastline  and  harbor  perimeter.  The  Test  1  grid  (36.106  elements  and 
18,483  nodes)  included  a  large  offshore  area  (da^ed  line  in  Figure  2)  and 
required  -10  min  of  central  processing  unit  (cpu)  on  the  Cray  Y-MP  for 
a  single  frequency  and  direction.  The  Test  2  grid  (6.911  elements  and 
3,679  nodes)  included  a  much  smaller  offshore  area  (dotted  line  in  Fig¬ 
ure  2)  and  required  only  -37  s  of  cpu  time  for  each  test. 

Comparisons  between  the  results  from  the  two  different  size  grids  (i.e.. 
Tests  1  and  2)  showed  only  small  differences  between  the  predicted  reso¬ 
nant  frequencies  and  sprttial  distribution  of  energy  at  those  frequencies. 
Although  the  resonant  frequencies  predicted  by  the  big  (Test  1)  and  small 
(Test  2)  grids  were  similar,  there  were  significant  differences  tetween  the 
magnitudes  of  the  resonant  peaks.  To  significantly  reduce  the  cpu  time  re¬ 
quited,  and  because  small  grid  results  are  adequate  for  determining  sensor 
placement,  subsequent  seiche  tests  for  this  study  primarily  used  the  small 
grid.  Future  work  comparing  the  magnitude  of  predicted  and  observed  res¬ 
onance  should  use  a  grid  large  enough  to  reduce  model  artifacts. 

Sea  surface  elevation  amplification  (with  respect  to  an  incident  wave 
of  amplitude  1.0  at  the  near-field  mode  boundary)  spectra  between  O.OOOS 
and  0.012  Hz  from  IS  locations  (solid  symbols  in  Figure  3)  within  the  har¬ 
bor  are  shown  in  Figure  5.  There  are  numerous  amplification  peaks  in 
this  frequency  band  with  varying  magnitudes  depending  on  the  location  in 
the  harbor. 

Examples  of  the  spatial  structure  at  several  lesorumt  fiequencies  are 
shown  in  Figures  6-9.  At  the  grave  mode  frequency  (-0.001  Hz,  Fig¬ 
ure  6).  relatively  uniform  amplifications  in  the  harbor  interior  decrease  to¬ 
wards  the  harbor  mouth.  This  very  low  frequency  mode  is  relatively 
broad  with  significant  amplifications  over  a  range  of  frequencies  (Fig¬ 
ure  S).  The  spatial  structure  of  tiw  narrower  peak  in  the  amplification 
spectrum  at  -0.003  Hz  (Figure  5)  is  shown  in  Figure  7.  This  seiche  rocks 
back  and  forth  between  the  east  and  west  sides  of  the  harbor  with  a  nodal 
line  (zero  crossing)  in  the  north-south  direction  in  the  middle  of  the 
harbor. 

At  higher  resonant  frequencies  the  spatial  patterns  of  wave  amplifica¬ 
tion  are  more  complex  with  multiple  nodes  and  antinodes.  Above 
0.006  Hz  the  amplification  spectrum  is  dominated  by  modes  with  edge- 
wave-like  spatial  patterns  (Figures  8-9).  These  modes  are  characterized 
by  nodes  and  antinodes  along  the  south  side  of  the  harbor  with  the  largest 
amplifications  along  the  shoreline  and  decreasing  amplifications  towards 
the  deeper  water  in  the  center  of  tlm  harbor. 
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Figure  5.  HARBD  model-predicted  sea  surface  elevation  amplification 
spectra  with  respect  to  incident  waves  of  amplitude  1.0  at 
the  near-field  boundary  (dotted  line  in  Figure  2)  for  inviscid 
waves  (Test  2)  at  various  locations  in  the  harbor  (shown  by 
solid  symbols  in  Figure  3) 


Constant  Depth  Harbor 

Model  tests  with  a  constant  depth  harbor  were  run  to  determine 
whether  the  ’edge’  modes  in  Figures  8  and  9  were  due  to  the  overall  shape 
of  the  harbor  or  a  result  of  refractive  trapping  of  energy  by  the  shallow 
depths  along  the  west  and  south  sides  of  the  harbor  (i.e.,  edge  waves).  In 
Test  3  the  offshore  depths  were  not  altered  but  the  harbor  depth  was  uni¬ 
formly  13  m. 

Amplification  spectra  in  this  case  (Rgure  10)  have  many  fewer  reso¬ 
nant  modes  than  with  the  field  bathymetry  (compare  Figures  5  and  10). 
None  of  the  resonant  modes  for  this  constant-depth  harbor  have  the  ‘edge’ 
spatial  patterns  predicted  with  the  real  bathymetry.  These  results  indicate 
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Figure  6.  Contours  of  model-predicted  wave  amplification  with  respect  to 
unit  wave  amplitude  at  the  near-field  boundary  (dotted  line  in 
Figure  2)  for  inviscid  waves  at  ^*0.0011  Hz  (Test  2).  Dashed 
lines  are  amplification  of  -0.  Circles  indicate  the  selected 
harbor  sensor  sites 

that  the  predicated  ‘edge’  modes  are  predominantly  a  function  of  the  shal¬ 
low  bathymetry  and  not  the  shape  of  the  harbor.  If  these  ‘edge’  modes  are 
energetic  in  the  Held  observations,  it  is  likely  they  could  be  suppressed  by 
dredging. 


Bottom  Friction 

Sensitivity  of  the  resonant  frequencies  and  spatial  patterns  of  amplifica¬ 
tion  to  bottom  friction  was  tested  using  a  constant  bottom  friction  coeffi¬ 
cient  h=0.001  (Test  4).  The  results  showed  only  slight  shifts  in  the 
resonant  frequencies  with  little  change  in  the  spatial  patterns  of  amplifica¬ 
tion  compared  to  the  inviscid  results  (Test  2).  However,  with  the  excep¬ 
tion  of  the  grave  mode  (-0.001  Hz)  which  showed  only  a  slight  reduction 
in  ampliHcation,  magnitudes  of  the  resonant  peaks  were  drastically  re¬ 
duced  by  bottom  friction. 
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Figure  7.  Contours  of  wave  amplification  for  f-0.0029  Hz  (see  caption  for 
Figure  6) 


Boundary  Reflection 

Sensitivity  of  the  harbor  amplifications  to  boundary  absorption  was 
also  tested  (Test  5).  The  results  showed  that  a  small  change  in  boundary 
reflectivity  (R=0.9S)  results  in  significant  decreases  in  the  amplifications 
but  has  little  effect  on  the  frequencies  and  spatial  structure  of  the  resonant 
modes. 


Rubble-Mound  Breakwaters 


The  purpose  of  the  last  test  (Test  6)  was  to  determine  the  sensitivity 
of  the  seiche  modes  to  possible  transmission  of  waves  over  the  rubble- 
mound  breakwaters.  Replacing  all  or  part  of  the  breakwater  with  a  sub¬ 
merged  breakwater  is  one  way  to  simulate  wave  transmission  through 
a  boundary  in  HARBD  (Houston  1976).  The  end  portions  of  the  break¬ 
water  (an  ~200-m-long  section  of  the  west  breakwater  and  an  --4()0-m-long 
section  of  the  east  breakwater)  are  relatively  thin,  rubble-mound  sections 
that  may  be  susceptible  to  transmission.  These  portions  were  removed  in 
the  mo^l  test  grid  by  replacing  the  zero-depth  breakwater  surface  area 
with  2-m  water  depths.  Although  the  reflecting  elements  of  the  breakwater 
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Figure  8.  Contours  of  wave  amplification  for  fsO.0065  Hz  (see  caption  for 
Figure  6) 

boundary  were  no  longer  present  in  the  model,  the  effect  of  the  break¬ 
waters  was  not  completely  removed  because  of  the  shallow  depths  which 
remained  at  the  breidtwater  location. 

Resonant  frequencies  and  spatial  patterns  of  amplification  for  this  har¬ 
bor  with  shallow  shoals  instead  of  breakwaters  were  similar  to  the  results 
obtained  for  the  inviscid  harbor  in  Test  2.  There  was  a  slight  shift  in  the 
grave  mode  peak  frequency  (-0013  Hz)  for  this  harbor  geometry  com¬ 
pared  to  the  harbor  with  impermeable  breakwaters  (i.e.  Test  2),  and  magni¬ 
tudes  of  amplification  were  also  different.  Some  modes  in  the  harbor 
without  bre^waters  are  more  resonant  (i.e.,  bigger  amplifications)  and 
some  are  less  resonant  than  the  harbor  with  breakwaters). 
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Figure  9.  Contours  of  wave  amplification  for  f>0.0077  Hz  (see  caption  for 
Figure  6) 
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Amplification 


Figure  1 0.  Model-predicted  amplification  spectra  (see  caption  for  Fig¬ 
ure  5)  for  a  constant  depth  {13-m)  harbor  (Test  3)  at  various 
locations  in  the  harbor  (shown  by  solid  symbols  in  Figure  3) 
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5  Site  Selection 


Harbor  Sensors 

Sites  chosen  for  the  four  harbor  sensors  are  shown  by  lettered  circles  in 
Figure  3  and  circles  in  Figures  5-8.  Amplification  spectra  from  HARBD 
Test  2  for  these  four  sites  are  shown  in  Figure  1 1 .  Sensor  A  will  be  be¬ 
tween  three  main  piers  where  most  of  the  present  harbor  operations  are 
located  and  sensor  B  will  be  positioned  in  the  vicinity  of  a  proposed  pas¬ 
senger  ship  and  ferry  terminal  (HDOT  1989).  Sites  A  and  B  are  expected 
to  be  dominated  by  the  ’normal*  harbor  seiche  modes  (Figures  6  and  7). 

Sensors  C  and  D,  located  along  the  shallow  southeast  perimeter  of  the 
harbor,  will  be  used  to  determine  the  signiHcance  of  the  predicted  ’edge’ 
modes  (e.g..  Figures  8  and  9).  Future  modifications  to  the  harbor  are  also 
being  proposed  for  the  area  in  the  vicinity  of  these  sensors,  adding  to  the 
importance  of  resolving  these  potentially  energetic  ‘edge’  modes.  For  ex¬ 
ample,  if  a  small  area  around  site  D  were  dredged  to  utilize  the  existing 
pier,  ‘edge’  modes  along  the  south  boundary  could  be  a  very  serious 
problem. 

Hme  series  of  2.3-br  duration  will  be  collected  at  harbor  sensors  A,  B, 
C,  and  D,  with  a  sample  rate  of  1 .0  Hz  (2.0  Hz  for  site  D),  which  will 
allow  for  the  resolution  of  sea,  swell,  and  infragravity  waves  at  these  sites. 


Offshore  Array  Site 

The  site  for  the  offshore  slope  array  was  selected  to  measure  both  en¬ 
ergy  levels  of  the  infragravity  frequency  motions  outside  the  harbor  and 
incident  sea  and  swell  representative  of  the  area.  The  selected  location 
was  highly  constrained  by  logistical  considerations.  Shipping  lanes  must 
be  avoided  for  obvious  reasons  and  sites  deeper  than  ~1 8  m  were  not  con- 
sittered  both  because  of  the  limited  bottom  time  for  scuba  divers  during  in¬ 
stallation  and  the  attenuation  of  high-frequency  seas  at  the  seafloor. 

Cable  lengths  also  limited  the  potential  array  sites. 
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Figure  11 .  Amplification  spectra  (see  caption  for  Figure  5)  for  inviscid 
waves  (Test  2)  at  the  selected  harbor  sensors  (shown  by  let¬ 
tered  circles  in  Figure  3) 


A  refraction-diffraction  numerical  model  (O’Reilly  and  Guza  1993) 
was  used  to  simulate  the  transformation  of  sea  and  swell  waves  from  deep 
water  to  offshore  of  the  harbor.  Examples  of  the  swell  (peak  period  =  17  s) 
height  (normalized  by  the  deepwater  height  Hq  offshore  of  the  harbor)  for 
two  incident  wave  directions  are  shown  in  Figure  12.  Model  results  indi¬ 
cate  that  the  offshore  bathymetry  at  this  site  focuses  the  incoming  waves 
resulting  in  convergence  zones  of  high  energy  at  certain  locations  which 
depend  on  the  incident  wave  direction.  Wave  energy  measured  in  these 
convergent  zones  will  be  site-speciHc  and  only  qualitatively  representa¬ 
tive  of  the  wave  conditions  over  a  larger  area.  None  of  the  potential  sites 
considered  was  out  of  these  convergence  zones  for  all  probable  incident 
wave  directions. 

The  site  selected  for  the  slope  array  is  shown  in  Figures  2  and  12.  At 
this  location  the  slope  array  will  at  least  n^asure  the  sea  and  swell  con¬ 
ditions  that  are  impinging  directly  on  the  harbor.  Results  from  HARBD 
using  the  large  offshore  grid  (Test  1)  also  indicate  that  this  slope  array  site 
is  not  located  near  nodes  of  the  resonant  haibor  modes.  Hme  series  of 
2.3-hr  duration  sampled  at  1 .0  Hz  will  be  collected  at  each  of  the  four 
pressure  gauges  in  the  slope  array. 
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Figure  12.  Refraction-diffraction  model  predicted  swell  (peak  period 

■  17  8)  heights  H  (normalized  by  the  deepwater  wave  height 
Hq,  scale  is  on  the  right-hand  side)  offshore  of  Kahului  Har¬ 
bor  for  incident  waves  from  the  north  (e-355*,  Figure  12a) 
and  northeast  (6a15*,  Figure  12b)  directions.  The  square 
symbol  is  the  site  selected  for  the  array.  Depth  contours 
shown  are  in  meters 
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Optional  Autonomous  Sensor  Site 


An  optional  fifth  nondirectional  autonomous  gauge  is  available  for  the 
Kahului  study.  The  gauge  is  an  internal  recording,  battery-operated  gauge 
that  could  be  placed  without  regard  to  cable  restrictions.  It  will  operate 
approximately  6  months,  depending  upon  the  sampling  scheme.  Four  pos¬ 
sible  regions  were  considered  for  this  gauge: 

a.  Inside  the  harbor,  to  better  define  harbor  oscillations  and  provide 
additional  calibration  for  the  harbor  models.  The  southwest  comer 
is  the  most  likely  site  within  the  harbor. 

b.  North  of  the  harbor,  to  better  define  the  incident  wave  field  and 
provide  additional  calibration  for  the  wave  transformation  model.  A 
site  on  the  western  side  of  the  entrance  channel,  to  better  define  the 
shadowing  effects  of  the  reef  on  the  eastern  side,  is  the  most  likely 
site. 

c.  East  of  the  harbor,  to  better  deflne  the  incident  waves  for  potential 
reclamation  on  the  shallow  region  seaward  of  pier  1 .  Data  from  this 
site  would  improve  selection  of  design  conditions  for  revetments. 

d.  Inside  the  western  breakwater,  to  better  define  the  transmission 
characteristics  of  this  stracture.  There  is  no  reliable  way,  using 
nondirectional  gauges  in  the  hartmr,  to  separate  energy  coming 
through  the  entrance  from  energy  coming  through/over  the 
breakwater.  Accurate  representation  of  the  permeability  of  these 
structures  may  be  possible  in  the  physical  models,  with  appropriate 
calibration.  A  gauge  inside  the  breakwater,  located  above  the 
normal  high-water  line,  will  record  the  effects  of  large  waves 
overtopping  the  structure. 

The  four  alternatives  are  listed  in  order  of  decreasing  confidence  in  the 
ability  to  effectively  model  the  processes  being  measured,  and  thus  make 
direct  application  of  the  measurements.  This  does  not  imply  they  are  in 
decreasing  order  of  value.  A  fifth  gauge  in  the  harbor  provides  incremen¬ 
tal  information  over  four  gauges,  but  given  the  complex  bathymetry  off¬ 
shore,  the  R/D  wave  transformation  model  would  benefit  even  more  from 
a  second  calibration  point.  However,  a  directional  gauge  would  be  more 
useful  for  thi.s  purpose  than  a  nondirectional  gauge.  Design  wave  informa¬ 
tion  is  needed  for  future  revetment  design,  but  unless  an  extreme  event  oc¬ 
curs  while  the  gauge  is  in  place,  it  will  be  difficult  to  predict  conditions 
on  the  reef  flat  from  statistically  derived  extreme  conditions  offshore.  Fi¬ 
nally,  the  need  to  accurately  model  the  processes  in  the  harbor,  which  in¬ 
clude  wave  overtopping  during  high  wave  conditions,  is  countered  by 
uncertainty  in  the  ability  to  reduce  the  complex  signal  from  a  gauge  in  this 
environment  to  a  usable  calibration  product  for  any  model. 
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Selectioa  the  best  locadoa  will  not  be  made  at  this  time.  Analysis  of 
the  data  from  the  existing  gauges  and  experience  with  die  models  will  pro¬ 
vide  the  best  guidance  on  the  need  for  actional  measurements.  If  indicated, 
the  optimal  site  for  an  additional  gauge  will  be  decided  later  in  the  study. 


Summary 

The  immediate  purpose  of  the  sensor  placement  study  was  to  design  a 
wave  measurement  program  which  would  aid  in  understanding  forcing  and 
damping  of  seiche  in  the  existing  harbor,  and  also  to  empirically  deter¬ 
mine  certain  harbor-speciHc  modeling  parameters  vital  for  the  study  of  fu¬ 
ture  modincations.  Limitations  on  resources  and  logistics  constrained  the 
experiment  design.  The  resources  allocated  to  the  wave  measurement 
study  include  a  directional  slope  array  to  measure  the  incident  wave  condi¬ 
tions  offshore  of  the  harbor  and  four  pressure  sensors  to  measure  the  har¬ 
bor  seiche.  An  additional  autonomous  (self-contained)  pressure  sensor  is 
available  and  possible  locations  for  this  gauge  were  discussed.  The  re¬ 
sponse  of  Kahului  Harbor  to  infragravity  frequency  (0.001  -  0.01  Hz) 
waves  was  numerically  modelled  with  tte  HARBD  model. 

Results  of  the  numerical  model  tests  show  that  Kahului  Harbor  is  char¬ 
acterized  by  numerous  resonant  modes  in  the  frequency  band  from  0.001 
to  0.012  Hz.  These  modes  can  be  qualitatively  separated  into  what  we 
refer  to  as  'normal*  seiche  modes  and  'edge*  modes.  The  'normal*  modes 
(e.g..  Figures  6  and  7)  are  a  function  primarily  of  the  hariior  shape  and  av¬ 
erage  depth  with  few  nodes  (zero  amplitudes)  within  the  harbor.  In  con¬ 
trast,  the  'edge*  modes  (e.g..  Figures  8  and  9)  are  a  function  primarily  of 
the  shallow  sloping  south  harbor  boundary  and  are  manifested  by  nodes 
and  antinodes  Signed  along  this  boundary  and  high  ampliHcations  along 
the  harbor  shoreline.  The  spatial  structure  and  frequencies  of  seiche  were 
not  sensitive  to  reasonable  variations  in  bottom  friction,  boundary  reflec¬ 
tion,  and  breakwater  transmissivity  as  prescribed  in  the  HARBD  numeri¬ 
cal  model. 

An  R/D  wave  transformation  numerical  model  was  used  to  estimate  off¬ 
shore  sea  and  swell  conditions.  Based  on  the  results  of  the  (uncalibrated) 
models,  sensor  placements  were  selected  which  are  best  suited  to  charac¬ 
terize  the  frequencies,  spatial  structure,  and  magnitude  of  the  resonant 
modes  and  which  quantify  the  most  important  modeling  unknowns:  (a)  the 
energy  level  of  low-frequency  motions  outside  the  harbor  and  their  vari¬ 
ous  driving  forces  (e.g.,  energy  of  sea  and  swell)  and  (b)  the  damping  (in¬ 
side  the  harbor  and  breakwaters)  which  can  dramatically  reduce  the 
amplitude  of  resonant  seiches. 

Recommendations  from  this  report  were  followed  in  deployment  of  the 
Kahului  Harbor  gauge  networb  in  October  of  1993.  The  network  was  in¬ 
stalled  and  is  operated  by  the  (3cean  Engineering  Research  Group  at  the 
SIO. 
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